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Abstract 

The paper is devoted to a systematic study of the duahty of processes in the 
sense that Ef{Xf,y) = Ef{x,Y^) for a certain /. This classical topic has well 
known applications in interacting particles, intertwining, superprocesses, stochastic 
monotonicity, exit - entrance laws, ruin probabilities in finances, etc. Aiming mostly 
at the case of / depending on the difference of its arguments, we shall give a sys- 
tematic study of duality via the analysis of the generators of dual Markov processes 
leading to various results and insights. 
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1 Introduction 



The paper is devoted to a systematic study of the duality of processes in the sense that 

Ef{Xf,y) = Ef{x,Y^) for a certain /. This classical topic has well known applications 

"-j ' in (and deep links with) interacting particles (see e.g. [25] and references therein), inter- 

rS ■ twining (see e.g. [Z], [H], [l3], [22]), superprocesses (see [I5], |26]), stochastic monotonicity 

_ci_- (see e.g. [32] and [10]), exit - entrance laws (see [H]), ruin probabilities in finances (see 

[12]). birth and death processes (see [31], [2]) and others. Aiming mostly at the case of / 

depending on the difference of its arguments, we shall give a systematic study of duality 

via the analysis of the generators of dual Markov processes (extending the analysis of 

one- dimensional processes from [22], [21]) leading to various results and insights. 

1.1 Objectives 

In stochastic analysis one meets various kinds of duality. For instance, the Markov pro- 
cesses Xf and Y^ with values in a Borel space X are called dual with respect to the 
reference measure u on X, if the duality equation 

f Eh{X^)g{x)u{dx) = f h{x)Eg{Y^)u{dx) (1) 

Jx J X 

holds for an appropriate class of functions h,g, see e.g. [3] and references therein for. 



In another approach, the Markov processes Xf and Y^ (small x, y here and in what 
follows stand for the initial points) with values in possibly different Borel spaces X and 
Y are called dual with respect to a Borel function / on X x F, if 

E/(Xf,y) = E/(x,r/) (2) 

for all X e X, y G y , where E on the left hand side and the right hand side correspond to 
the distributions of processes Xf and Y^ respectively, see e.g. [25] and references therein 
for an extensive application of this notion in interacting particles. 

A particular case of (|2]) is the duality of one-dimensional processes (X and Y are 
real-valued) arising from stochastic monotonicity, where f{x,y) = l{x>j/} (we denote here 
and in what follows by lj\/ the indicator function of the set M) and hence ([2]) turns to 
the equation 

P(Xf >y) = P(F/ < x), (3) 

see |29]. Other useful cases include f{x,y) = e^^ or f{x,y) = x^, used in particular in 
the theory of superprocesses, see e.g. Ch. 4 of [15] or Ch. 1 of [H]. For an application of 
duality in actuarial science see [T^ . 

The analytic analogs of the duality of the 1st kind is successfully used in the theory of 
operator semigroups independently of their probabilistic content, see e.g [3] and references 
therein. We shall start now with a sketch of a systematic study of duality obtained by 
extending (|2]) to general purely analytic setting aiming at the extension of the theory of 
dH]) to dualities generated by partial orders and more general translation invariant dualities 
arising from / depending on the difference of their arguments. 

There are many applications of duality in population dynamics, branching processes 
and other areas, see e.g. pj and references therein. 

1.2 On the general notion of semigroup duality 

For a topological (e.g. metric) space X we denote by -B(X) and C(X) the spaces of 
bounded Borel measurable and bounded continuous functions respectively. Equipped 
with the sup-norm ||/|| = sup^, |/(x)| both these spaces become Banach spaces. Bounded 
signed measures on X are defined as bounded a-additive functions on the Borel subsets 
of X. The set of such measures Ai{X) equipped with the total variation norm is also a 
Banach space. The standard duality between -B(X) and Ai{X) is given by the integration: 

(/,/^) = / f{x)fx{dx). 



X 



Let X, Y be two topological spaces. By a signed (stochastic) kernel from X to V 
we mean a function of two variables p{x,A), where x G X and A are Borel subsets of 
Y such that p{x, .) is a bounded signed measure on Y for any x and p{.,A) is a Borel 
function for any Borel set A. We say that this kernel is bounded if sup^ \\p{x, .)|| < oo. 
We say that this kernel is weakly continuous if the mapping x i— ?■ p{x, .) is continuous with 
measures Ai(Y) considered in their weak topology. If all measures p(x, .) are positive, 
the corresponding kernel is called a stochastic kernel. 

Any bounded kernel specifies a bounded linear operator B{Y) — ?■ B{X) via the formula 

Tg{x) = / g{z)p{x,dz). 



T is said to be the integral operator with the kernel p. The standard dual operator T' is 
defined as the operator M{X) — )• J\4(Y) specified by the duality relation 

(/,T'/i) = (T/,/i), 

or explicitly as 



T' ii[dy) = I p{x , dy) fi{dx) . 
Jx 

Clearly the kernel p{x, dz) is weakly continuous if and only if T acts on continuous 
functions, that is, T : c\y) -^ C{X). 

Let p{x, dz) be a bounded signed kernel from X to itself, T the corresponding integral 
operator, and let f{x,y) be a bounded measurable function on X xY. Let us say that 
the operator T^^^) : fi(r) -^ B{T) is f-dual to T, if 

{T''^^^f{x,.)){y) = {Tf{.,y)){x) (4) 

for any x, y, that is, the apphcation of T^ to the second argument of / is equivalent to 
the application of T to its first argument. Of course, if T^^^^ is /-dual to T, then T is 
/-dual to T^^^^ with f{y,x) = f{x,y). 

We say that / separates points of X if, for any xi, X2 G X, there exists y eY such that 
/(xi, y) 7^ /(a;2, y). The following is a bit more nontrivial notion. We say that / separates 
measures on X if, for any Qi, Q2 G M.{X), there exists y G y such that J /(x, y)Qi{dx) 7^ 
/ f{x,y)Q2{dx). If this is the case, the integral operator F = Ff : Ai{X) — )■ -B(F) given 
by 

{FQ){y) = j f{x,y)Q{dx) (5) 

is an injective bounded operator, so that the linear inverse F~^ is defined on the image 
F{M.{X)). Let us say that the function FQ is /-generated by Q. 

Remark 1. In /I7]/ the authors call a function g to be representable by f , if there exists 
a unique Q such that g = FQ. Paper /Ti] / deals with the application of duality to exit and 
entrance laws of Markov processes. 

1.3 Basic tools 

Proposition 1.1. Let f be a bounded measurable function separating measures on X and 
T an integral operator in B(X) with a bounded signed kernel p. Then T'^^^' is well defined 
on F{M.{X)) and its action on the f -generated functions coincides with T' , that is 

T^U) =FoT' oF-\ (6) 

or equivalently 

^-1 o T^(/) = r o F^^ (7) 

In other words, the f-dual operator T^^^^ is obtained by the 'dressing' of the standard dual 
T' by the operator F . 



Proof. Let g G F{M.{X)) be given by g{y) = f f{x,y)Qg{dx). Then 



X 



{Tf{.,y)){x)Qg{dx) = / / f{z,y)p{x,dz)Qg{dx) = / f{z,y)Q{dz), 
X Jx Jy Jy 



with 



Q{dz) = p(x,dz)Qg{dx). 

Thus T^^-f^g is /-generated hy Q = T'Qg, as required. D 

Remark 2. Equation dTj) zs a particular case of the so-called intertwining, see j^, /^, as 
well as fTSf . f^ , fT8f for exciting recent developments. Relations ([7]) for discrete Markov 
chains are analyzed in detail in U9^ . 

Representation ([6]) has a direct imphcation for the theory of semigroups. 

Proposition 1.2. Let f be a bounded measurable function separating measures on X and 
Tt a semigroup of integral operators in B{X) specified by the family of bounded signed 
kernel pt{x, dz) from X to X , so that Tq is the identity operator and TtTg = Tt+s, which, 
in terms of kernels, rewrites as the Chapman- Kolmogorov equation 



Pt{x, dz)ps{z, dw) = pt+s{x, dw). 
' x 

Then the dual operators T^ in F(Ai{X)) also form a semigroup, so that 

T,^(^) = F o t; o F-^ (8) 

Proof. This is straightforward from ([6]) and the standard obvious fact that T/ form a 
semigroup in Ai{X). D 

Remark 3. The duality ([T]) is of course also included in the general scheme above, that 
is, the dual can still be expressed as (El). For instance, if v{dx) has a density v{x) with 
respect to Lebesgue measure and T' can be reduced to the action on functions, then F~^ 
is the multiplication on nu{x) and f{x,y) = 6{x — y)v~^{x). 

It is also worth noting that the assumption of boundedness of / is not very essentiaL If 
it is not bounded (and we shall discuss interesting examples of such situations later), the 
integral operator F will not be defined on all bounded measures, but only on its subspace. 
This will be reflected in the domain of T^^^^ , but the whole scheme of Proposition 11.11 
still remains valid. 

1.4 Links with differential equations and stochastic processes 

Let us explain briefly the main ideas on the application of the above results to the theory 
of differential equations and stochastic processes. Precise details for concrete situations 
will be discussed below. 



Let a semigroup Tt in B{X) be generated by a (possibly unbounded) operator L in 
B{X) defined on an invariant (under all Tt) domain D C B{X), so that 



d 
It 



Tth = \im-(Tth - h) = Lh, he D, 



with convergence in some appropriate topology (say, strongly or point-wise) and thus the 
operators Tt represent resolving operators for the Cauchy problem of the equation h = Lh. 
Then ([6]) implies that 



d_ 
di 



1. ^ 'q = h o —- 
df 



r o F-^g = FoL'o F-^g, 
t=o 



that is, the generator of the semigroup T^ is 

L^(/) =FoL'oF-\ (9) 

so that the operators T^ represent resolving operators for the Cauchy problem of the 
equation g = L^^^^g. Here L' is of course the standard dual operator to L. Thus duality 
can yield explicit solutions for equations of this kind. Of course, our arguments were 
heuristic as we did not pay attention to the domain of definition of L', which should be 
done in concrete situations. The main difficulty here is to characterize the operator Ff. 

Next, in order to be able to fill the duality equation @ with probabilistic content, 
i.e. to rewrite it as 02]), the semigroups Tt and T^ should be positivity preserving and 
generate some Markov processes. 

This question effectively reduces to the question of whether, for a given conditionally 
positive operator L, the corresponding dual L^'^^^ is also conditionally positive. 

It is seen now that the basic questions to be addressed to make the theory work for 
concrete functions / are (i) the characterization of the operators F and F~^ (for the 
analytic part of the story) and (ii) the criteria for conditional positivity of L^*^-^-* (for its 
probabilistic content). 

As we shall see it is often convenient to reduce the operator F to some subclass of 
Borel measures Q, where its inverse can be explicitly found. For instance, it is often easier 
to work with Q having density with respect to some reference measure. 

1.5 Content and plan of the paper 

We shall apply formulas ([H]) and (jH]) to characterize classes of dual Markov processes with 
respect to various functions / depending on the difference of its arguments. Section [2] deals 
with duality on R*^ arising from Pareto and similar partial orders. The full characterization 
of duality is given in terms of generators for basic classes of Feller processes. Section [3] 
discusses several examples of duality with operator F~^ being the Laplacian or a fractional 
Lapacian. Section H] initiates an application of formulas (jH]) and (jH]) to the study of duality 
for processes in domains with a boundary. To circumvent specific difficulties arising from 
the boundary, we introduce here an additional tool of a regularized dual. 

The extension of the theory to time-nonhomogeneous Markov processes will be ana- 
lyzed in 



2 Duality from orders and other binary relations 

2.1 Basic notions 

As our basic example we consider /-duality for functions / arising from translation- 
invariant partial orders, or more generally, from translation-invariant binary relations. 
Namely, let X be a topological linear space and M a Borel subset of X. Then M defines 
a translation-invariant binary relation Rm on X such that xRmV means, by definition, 
that X — y E M, or x E y + M. 

Let M = {{x,y) E X x X : xRMy}- Let us say that the duality (j4]) arises from the 
binary relation M, if 

fix, y) = fnix, y) = iMix, y) = Ix-niy) = ly+Mix). (10) 

Remark 4. // f -duality arises from a translation-invariant binary relation Rm and if 
both Tt and T^ are known to be integral operators with kernels pt{x, dz) andp^ {y, dw) 
respectively, one can give another instructive proof of Proposition \l.<l\ bypassing represen- 
tation ([6]) and using instead Fubbini's theorem, as was done in 129^ for standard one- 
dimensional duality. Namely, it is sufficient to show the semigroup identity T^+/ = 
Ts Tf- applied to the functions /(x, .) = 1x-m, o,s it then extends to the whole 
F(A^(X)) by linearity. And for these functions we have 

{T^+Plx^M){y) = (Tt+slj/+M)(x) = {Tt{Tsly+M))ix) = / pt{x,dz){Tsly+M)iz) 



Ptix,dz)iTf^^h,.M)iy) = j Pt{x,dz) U l,.M{w)pf^f\y,dw)\ . 
Applying Fubbini's theorem this rewrites as 

pf^f\y,dw) [ 1^+M{z)pt{x,dz) = [ {Ta^+^^{x)p^^f\y,dw) = T^^f\T['^fh,.M){y), 



as required. 

If M contains the origin and is closed under the addition of vectors, then the relation 
Rm is a pre-order (i.e. it is reflexive and transient) and can be naturally denoted by >m- 
If this is the case and Tt and Tj are integral operators with positive stochastic kernels 
thus specifying Markov processes, then duality relation (jl]) or equivalently (E]) turns to 
the equation 

V{X^>My) = V{Yy<Mx), (11) 

extending one-dimensional duality ([3]). 

The basic example we are going to analyze now is the Pareto partial order in X = R"^, 
i.e. >M with M = R!^, and its natural extension with M = C(ei, ■ ■ ■ , e^) the cone 
generated by d linear independent vectors {ei, ■ ■ ■ , e^} in R'^: 

d 

C{ei,--- ,ed) = {x = ^ajej : a^ > 0, j = 1, ■ ■ ■ , c?}. (12) 



Of course the relation >m with such M is again a Pareto order, but in a transformed 
system of coordinates. 

Let us start with M = R^ corresponding to the Pareto order, which we shall denote 
just by > omitting the subscript M. The corresponding dual semigroups or processes (if 
exist) will be referred to as Pareto dual. In this case 



iFQ)iy) = / fMix,y)Qidx) = / Qidx) (13) 

is just the usual multidimensional distribution function for the measure Q on R"^. It is 
known (and easy to see) that FQ characterizes Q uniquely implying that F is injective 
and thus /m separates measures on R'^ yielding the main condition of Proposition 11.11 
Moreover, if Q has a density q with respect to Lebesgue measure, then q can be found 
from FQ = g hy differentiation: 

q{yi, ■■■,y,) = F^'giy) = (-1)'^^^^^. (14) 

Thus, for the Pareto order, the operator F^^ has the simple explicit expression. 

In the case of the orders arising from the cones M = C{ei, ■ ■ ■ , e^) given by (IT2l) this 
formula generalizes to 

q{yi, ■■■ ,ya) = {F 'g){y) = -I^ttt ^, 15 

|det(ei,e2,- ■■ ,ed)\ 

where det(ei,e2, • ■ ■ , e^) = det(e]) is the determinant of the matrix whose ith columns 
consist of the coordinates of the vector e,- and 



d'^'g sr-^ d'^g 



o*2 „«d 



Remark 5. For completeness, let us sketch a proof of this formula. If a measure Q on 
R*^ has a continuous density q, so that 



g{x) = FQ{x) = / q{z)dz, 

ly+C{ei,-,ea) 



the function q can be clearly found as the limit 

q{y) = lim q{z) dz\U{hei,- ■ ■ ,hed)\~^, (16) 

'^~^^ Jv+n(hei.---.he^) 



'y+n(hei,--- ,hed) 

where 



n(/iei, ■ ■ ■ , hcd) = {x = y^ ctjhcj, aj E [0, 1]} 



is the parallelepiped built on the vectors {hci, ■ ■ ■ , hcd} and 

\U{he,,--- ,hed)\ = h''\det{el)\ 
is its Euclidean volume. 
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From simple com,binatorics it follows (see e.g. 12^ ) that 

q{z) dz 

y+Tl{hei,--- ,he^) 

= 9{y) -^g{y + he,) + ^g{y + het + he^) + ■■■ + {-ifgijj + /id + ■ ■ • + hed). 

Let us expand all term,s in Taylor series up to the derivatives of order d. As the final 
expression should he of order h'^ (to get a limit in ( TT6|) ) we conclude that all terms with 
the derivatives of orders less than d necessarily cancel, so that 

I q{z) dz 

Jy+U{hei,-- ,heci) 

=>'(-i:0N+E0i^.+«j+-+(-i)'0i«.+--+'j)+o(ft«). 

(17) 
where 0{h'^^^) denotes the expression of order h"^^^ that does not contribute to the limit 
in (iT6l) . and where we use the well established (though a bit ambiguous) notation for the 
action of the higher order derivative on equal vectors: 

It remains to note that all terms in expansion (fTTl) containing products of coordinates of 
coinciding vectors should vanish (otherwise, using different scaling on Cj we would arrive 
to a contradiction with the existence of the limit in (iT6l) ). The only non-vanishing terms 
should contain the products of d coordinates of all d vectors. All these products comes 
from the last term in the sum flTTl) leading to flTSl) . 

For instance, let us consider a 'two-dimensional light cone' 

C(el,e2) = {(a;,y):y>|x|}GR^ (18) 

corresponding to vectors ei = (1,1), 62 = (—1,1). Then formula (TT5l) for the inverse 
operator turns to the simple wave operator 

q{x, y) = F-^g{x, y) = - (^— - —J (x, y). (19) 

2.2 Duality from Pareto order: global analysis 

Let us now make the detailed analysis of the duality arising from the standard Pareto 
order in R'^, i.e. with M = R^. We aim at (i) finding explicitly the dual operator L^*^-^^ 
for the main classes of the generators of Feller processes in R*^ including diffusions and 
jump processes and (ii) establishing criteria (in terms of the initial operator L) ensuring 
that this dual operator is conditionally positive and specifies a Markov process, so that 
the duality relation ( TTT1) holds that we shall write simply as 

P{X^ >y)= P(r/ < x) (20) 



for the case of the Pareto partial order. 

Let us analyze formula ([6]) from Proposition II. 1[ In the case of duality arising from 
Pareto order and the operator T being integral with a probability kernel p{x, dz) (i.e. 
all measures p(x, .) are probability measures, as is the case for transition operators of 
Markov processes) it states that for a distribution function g oi a. measure Q on R'^. i.e. 

T^^i^g{x)=FoT' oF~^g{x)= f f p{z,dy)Q{dz). (21) 

We are interested in the question of when this operator can be extended to all bounded 
measurable g a.s a. positive operator preserving constants, i.e. as an integral operator with 
a probability kernel. 

Assume first that the measure Q has a continuous density q so that (TT4l) holds, i.e. 

q{x) = (-1)^ ^ 



dxi ■ ■ ■ dxd 
In this case 

T^^f^gix) = {-ir [ [ p{z,dy) ^^\ dz. (22) 

Jy>x JRd OZi ■ ■ ■ aZd 

We like to get rid of the derivatives of g. To be able to do it, let us assume that the 
kernel p(x, dz) is weakly continuous and has weakly continuous mixed derivatives, that 
is, for any / C {1, ■ ■ ■ , c?} (including {1, ■ ■ ■ , c?} itself) the mixed derivative 

^{z,dy) (23) 

OZi 

is a well defined weakly continuous kernel (possibly signed). Then, integrating the integral 
over z in (l22l) by parts d times and assuming that all boundary terms vanish, we get 

T''(/)«W^£(.W£^-*1^(...,)).. (24) 

This is an integral operator with the integral kernel (more precisely its density) 

f C/T) 

p^{x,z)= --{z,dy). 

Jy>x ozi--- dZd 

For this operator to be positive and constant preserving, necessary conditions are that, 
for all X G R'^, 

"^^ {z,dy)>0, (25) 

^^ {z, dy)] dz = 1. (26) 



ly>x 9zi ■ ■ -dZd / 

From the integration by parts it is seen that for the last condition to hold it is sufficient 
to assume that for any subset / C {1, • ■ ■ , (i} excluding the whole set {1, ■ ■ ■ , ci}, 

f f Qp\I\ 

lim / dzi I — — {zi,zj,dy) = 0, (27) 

Rl^l Jy>x C^/ 

9 



Zt—^ — OO 



and there exists a finite limit 

f f dp\^\ 

lim / dzi / — — {zj,zi,dy), (28) 

^/^^JrI^I Jy>x ozi 

which equals 1 for the empty set /. Moreover, one sees by inspection that this condition 
also ensures that integrating by parts (122]) for a g having finite density f ll4p . all boundary 
terms will in fact vanish, justifying equation f l2^ . 
Thus we have proved the following statement. 

Proposition 2.1. Suppose an integral operator T in i?(R'^) is given by a probability 
kernel p{x, dy) having all mixed derivatives (!23|) well defined and weakly continuous and 
such that (125 p holds, (12 7p holds for any subset I C {1, ■ ■ ■ ,d} excluding the whole set 
{1, ■ ■ ■ , d}, and there exists a finite limit fl2S]) . which equals 1 for the empty set I. Then 
the Pareto dual operator T^^-^^ is also an integral operator with a probability kernel. 

Condition (!25|) is of course not directly verifiable. Therefore we shall see how it can 
be read from the generator of the process. 

2.3 Duality from Pareto order: deterministic and diffusion pro- 
cesses 

We plan now to find the generators of the dual processes, when they exist. Let us start 
with the simplest case of deterministic processes generated by the first order differential 
operators of the form 

L0(x) = (&(x), V0(x)) = J2bA^)^- (29) 



In this case the dual operator is well defined on functions and 

dx 



d 

L'g{x) = -diy{gb){x) = -^^^[bjix)g{x)]. 



For a vector x = {xi, ■ ■ ■ , Xd) € R'^ let us denote Xi the vector in R''^^ obtained from 
X by deleting the coordinate Xi. For a function g{x) let us write g{zi,Xi) for the value of 
g on the vector, whose ith coordinate is Xi, and other coordinates are those of the vector 
z. Let us write dzj for the product of differentials dz^ with all fc = 1, ■ ■ ■ ,d excluding j. 

Integrating by parts and assuming that g decays quickly enough so that the boundary 
terms at infinity vanish, we have 



L^^f^gix) = FL'F-'g{x) = (-1)^+^ / E 7^ 



3 '' 



b.iz)- ^ '^'^ 



dzi ■■■dzd 

-^y^ j_ ^_ b^{Zj,Xj)^-^-^J,Z,,X,)dZj. 



dzi ■ ■ ■ dzd 



(30) 

In general one cannot simplify this expression much further, and this is not a condition- 
ally positive operator (it does not have a Levy-Khintchin form with variable coefficients) 
without further assumptions. 
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Proposition 2.2. Let L have form (!29|) with all bj G C^(R'^) (the space of bounded 
continuous functions with bounded continuous derivatives). Then L^^^^ is given by ( 130|) . 
so that the solution to the Cauchy problem of the equation g = L^^^^g is given by the 
corresponding formula (jS]) with F and F~^ given by fll3p and f ll4p . Moreover, if each bj 
depends only on the coordinate Xj, then 



LDU) 



9[X) 



'bj[Xj 



^ dg 
dxi 



(31) 



that is, L^^^^ coincides with L up to a sign and the dual process exists and is just the 
deterministic motion in the opposite direction to the original one. 

Proof. Formula ( 13T]) is straightforward from (|30l) and the assumptions made on bj. This 
makes the last statement plausible. However, strictly speaking, having the generator 
calculated on some subclass of functions does not directly imply that the semigroup T^^^^ 
coincides with the semigroups on C(R'^) generated by operator fl3Tl) . The simplest way to 
see that this is in fact the case is via durect calculations with the semigroup T^ itself. 
Namely, if the deterministic Markov process X^ with generator (!29|) can be expressed as 
X^ = X*(x) via the solutions X*(x) of the Cauchy problem for the ODE x = b{x), its 
transition kernel takes the form pt{z, dy) = 6{y — X^{z)). Then f l22p becomes 



T^'^'^gix) 



dg" 



-If / ^r^^dz. (32) 

Under the assumption that bi depend only on Xj, the coordinates of X^{z) are themselves 
solutions Xl{zi) of the one- dimensional ODE Xi = bi{xi), so that one has 



T, 



D{f) 



9[x) 



dg" 



-dz. 



xi(zi)>xi dzi- ■ ■ dzd 
From the obvious monotonicity of one-dimensional ODE this rewrites as 



rpD(f) 



9[x) 



' z,>{xiy^(x,) dzi---dzd 
which is of course the semigroup generated by the operator (13T11 . 
Let us turn to a diffusion operator having the form 

d 



dz = g{{r)-\x)), 



L(j){x) = {a{x)'V ,'V)(j){x) = 2_. (^ij{^) 



d' 



with a positive definite diffusion matrix a{x) 
In this case 



{a,j{x)). 



dxidxj 



{X 



(33) 

(34) 
D 

(35) 



d^ 



L'gix) = Y^ 7ir7^[%(^)^(^)]' 



dxidxj 



and consequently 

L^(/)^(a;) = FL'F-^g{x) 



■ly 



z>x 






d^ 



dzjdz 



J I 



aij[z) 



d'gjz) . 
dzi-- ■ dzd ' 



dz\ • • • dZd- 
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Let us integrate twice by parts the terms containing mixed derivatives and integrate once 
by parts the remaining terms. This yields 






^jA^J-i ^3 i 



dzi-- ■ dzd 



^Zj, Xj) 



dz-i 



d'^g 



dzi--- dzd 



\Zij , Xi , Xj J ClZij , 



where Zij denotes the vector in R obtained from z by deleting ith and jth coordinates, 
and {zij,Xi,Xj) is the vector with ith and jth coordinates taken from the vector x, and 
other coordinates taken from the vector z. In case d = 1, the second sum in this expression 
is of course empty. 

Again in general case one cannot simplify this expression essentially. However, as- 
suming additionally that the coefficients Ojj depends only on the coordinates Xi,Xj (in 
particular, an depends only on Xi), we have 



L^(f^g{ 



Xj 



d „ 



_d_ 

dxi 



+2( 



Kj 



(^ij \-^i ; -^j I 



ajj[Xjj 



Qdg 



dzi ■ ■ ■ dza 



\Zj, Xj j 



dzA 



dzi- ■■ dzd 



^Zij •, Xi^ "^ j } aZ{j. 



Integrating by parts with respect to the variables Zj in the first sum and the variables Zi. 
in the second, yields (assuming the boundary terms at infinity vanish) 



L^(/)^( 



Xj 



^ dxi 

7 = 1 ■' 



^jj \^j t 



dgjx) 
dxi 



~r ^ / ^ (^ijyXij Xj) 



i<j 



d'g 



or 



L^^^hix) = Lg{x) + J2 



dajj [Xj ^ 



dg{x) 



dxi 



dXj 



(36) 



(37) 



Proposition 2.3. Let L have form ( l35l) with a positive definite diffusion matrix a{x) = 
{aij{x)) and with all aij G C^(R'^), so that L generates a Feller diffusion in R*^ that we 
denote Xf . If the coefficients a^ depends only on the coordinates Xi,Xj, then L^^f^ is 
given by ( l37l) and it also generates a diffusion process in R'^ that we denote Y^ , and the 
duality relation ( l20l) holds. 



Proof. Again formula (l37j) makes the statement very plausible, but to deduce ([6]) from 
(|9]) additional argument is of course needed. This goes as follows. 

But notice first that it is sufficient to prove the statement under additional assumption 
that coefficients aij are infinitely smooth with all derivatives bounded (actually we need 
twice differentiability for the above calculation of L^^^"* and d times differentiability for the 
formulas of Proposition 12.11 to make sense) and the operator L is strictly elliptic, because 
any L of type (l35l) can be approximated by the sequence of L of the same form but strictly 
elliptic and with smooth coefficients. Passing to the limit in the duality equation allows 
one to prove its validity for the general case. 
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Next, under this smoothness and non- degeneracy assumption, it is well known from 
the standard theory of diffusions (or Ito's processes) that operator (1371) generates a unique 
Feller process such that its semigroup Tj preserves the space C^(R'^) of twice continu- 
ously differentiable functions vanishing at infinity with all its derivatives up to order two. 
Hence, the Cauchy problem for the equation 

9 = L^^'h 

is well posed in classical sense for initial functions go from C^(R'^). It is then straight- 
forward to see (j9]) that both functions T^ 5^0 and F o T/ o F^^qq satisfies this equation. 
Consequently these two functions coincide implying (|6]) for the semigroups Tj and T^ ^ \ 
as required. D 



Thus we have shown that under appropriate assumptions the /-dual operators to the 
first order and diffusion operators respectively are again first order and diffusion operators 
respectively defining the /-dual or Pareto dual processes. 

It is instructive to see which diffusions are self-dual. This is given by the following 
result that is a direct consequence of Propositions 12.31 and 12.21 

Proposition 2.4. Let 

'^ d'^d) 1 '^ d(J) 

L(t){x) = Y, ""^^-^^^'^J^^^^^^^) + 2 5Z^^i^^^)^(^) (38) 

ij=l J jf=l J 

with a positive definite (possibly not strictly) diffusion matrix a(x) = {aij{x)) such that 
ttij depend only on Xi,Xj and are continuously differentiable (with bounded derivatives). 
Then the diffusion generated by L is self-dual in the Pareto sense. 

2.4 Application to other cones 

Generalization of our results to orders arising from cones C(ei, ■ ■ ■ , e^) can be obtained 
by the change of variables, though the calculations quickly become rather cumbersome. 
Let us consider only the simple example of the two-dimensional cone fll8j) . The question 
we are going to answer is as follows: under what conditions the diffusion operator 

d'^g d^g d^g 

Lg{x,y) = a{x,y)—^ + 2b{x,y)——- + c{x,y)—^ (39) 

ax^ oxoy oy^ 

generates a diffusion that has a dual in the sense of the order generated by C, and how the 
dual generator looks like. Having in mind the relation with the standard Pareto order we 
can expect that the coefficients should depend in certain way on two arbitrary functions 
of one variable and one arbitrary function of two variables. This is in fact the case as the 
following result shows. 



Proposition 2.5. Let L of form fl39|) with smooth coefficients generate a Feller diffusion 
Xf . If the coefficients have the form 

a(x, y) = a{x + y) + (3{x - y) + uj{x, y), 

c(x, y) = a{x + y) + l3{x - y) - uj{x, y), (40) 

6(x, y) = a{x + y) - l3{x - y) 
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with some smooth functions a,P,u, then Xf has the dual diffusion Y^ so that ( ITT]) holds 
with M = 0(61,62) of form (ITSj) . where Y^ is generated by the operator 

L^^f)g = Lg + A{a'{x + y) + (5\x-y))^{x,y) + A{a\x + y)-(5\x-y))^{x,y). (41) 

Proof. Formulas ( HOl) are obtained from Proposition 12.31 by rotation of coordinates, that 
is by change x' = x + y,y' = x — y. D 



2.5 Duality from Pareto order: jump processes 

Let us now turn to the generators L of pure jump processes, that is 



L4>{x) 



R'* 



{(p{w) — (t>{x))h'{x, dw) 



(42) 



with some bounded stochastic kernel z/. For a measure Q having a density with respect 
to Lebesgue measure, let us write shortly L'q for the measure L'Q. We have 



L'q{dz) = / q{x)v{x,dz)dx — q{z)dz / v{z,dw). 
Consequently, relabeling the variables of integration, we have 



FoL'{q) = {-lY f {L'q){dz) 

J z>y 



q{z)u{z, dw) dz. 



= (-1) / / q{z)u{z,dw)dz-{-l) 
The integrals in the two terms partially cancel. Namely, we can write 



F o L\q) = {-l)" I q{z) fl,>y f u{z, dw) - / u{z, dw) 

J \ \_J w>y 



+ l^^j/ / v{z,dw) ) dz, 

w>y 



implying 



F oL'{q) = {-l)'^ j q{z) l,^y u{z,dw)-l,>y u(z,dw) 

Hence, for a smooth {d times differentiable) function g we can write either 

d'^g{z 



dz. 



J OZi ■ ■ ■ OZd Jw>y Jw^y 



u{z, dw) 



dz. 



(43) 



or 



L^^f^g = FoL'oF~'g{y) 



(-l)M j ^^^u(z,dw)dz-{-l)' j [ -f^-uiz,dw)dz. (44) 
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If z/(z, dw) depends smoothly on z, this expression can be rewritten by moving the 
derivatives from g to u. For this transformation expression (H4|) is more handy than (H3|) . 
To perform the integration by parts in its second term we shall use the following simple 
formula (with a straightforward proof by mathematical induction) 



z>y 



dzi ■ ■ -dZf 



-(l){z)dz = {-if ^ / g{yj,zi)— — {yj,zi)dzi, 



(45) 



which is valid when the boundary terms at infinity vanish, for instance if either (f) or 
g vanish at infinity with all its derivatives. Here |/| is the number of indices in J, the 
integral over the set {zj > yj} is |/|-dimensional and {yj,zi) denotes the vector whose 
coordinates with indices from / are those of the vector z and other coordinates are from 
the vector y. 

Using this formula we transform (jH]) into the expression 



9[z, 



L''^f^g{y) = FoL'oF-^g{y) 
-—{z,dw)dz- y2 / dzigiVhZi) j — — {yj,zi,dw). 

w>y Jnd C^l ■ ■ ■ OZd /c|l~^ d} -^""l^yi -^^"^ -^ 

Singling out from the sum the terms corresponding to / being empty and I being the 
whole set {!,••• ,d}, this rewrites as 



d'^u 



9{z 



>y J^d dzi--- dzd 



[z, dw) dz 



d'^u 



9{z 



>y JRd dzi--- dZd 



{z, dw) dz 



~5Ilr n .1 / dzjg{yj,zj) / — — {yj, zj,dw) - g{y) / u{y,dw), 

where ^ denotes the sum over all proper subsets /, i.e. all subsets / excluding empty 
set and the whole set {1, ■ ■ ■ ,d}. Performing the cancelation in the first two terms yields 
finally (see the trick leading to fH5]) ) 

L^(^)(7(y) = FoL'o F-'g{y) = -g{y) [ u{y, dw) 



V / dzig{yj,zi) / — — 

^^ 'zi>yi JRd ozi 



/c{i,-,4 



+ / g{.z)dz 

'Rd 



'-ty 



d'^u 



iVh zi, dw) 



w>y 



dzi--- dzd 



[Z,dw)-U>y / — 

Jw^y ^Zi ■ ■ ■ 



dZrl 



[z, dw) 



(46) 



For instance, for d 



/v r f~i /*oo r ^ r 

g{z)dz —{z,dw)- g{z)dz —{z,dw)-g{y)u{y,dw), 

-oo J w>y ^^ J y J w<y ^^ J 

(47) 
which is the formula essentially obtained in [2T] and |22], and for d = 2 



L^^^^g{y) = -g{yi,y2) / v{y,dw) 
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-/ 9{zi,y2)dzi -—{zi,y2,dw)- g{yi,Z2)dz2 j ^—{yi,Z2, 

+ I ,(z„ Z2)dz,dz2 1.^, £^^ ^(., dn^) - 1.,, £^^ ^(., dw 



dw) 



{Ai 



Remark 6. It is worth stressing that one should he cautious in using these formulas as 
they may not he true for f not vanishing at infinity, say even for a constant function f (so 
that these formulas cannot he used even for checking conservativity condition L^^^^l = 0). 
Generally one has to use the following extension of (H5|) (also proved hy direct induction) 
that is valid whenever g, cf) are smooth and such that for all I G {!,■ ■ ■ , d} and yj there 
exist finite limits of the functions g{yi, zj), (f){yj, zj) and their derivatives in zj, as zj ^- oo 
(here oo means precisely +ooj.- 



z>y 



d''g{z) 
dzi ■ ■ ■ dz, 



-(f){z)dz 



-1/ E / 

IC{1,- ,d} '-y' 



d\I\, 



X](-l)'"^'^(l/JV, OOJ, Zl)-^ (l//\j, OOj, Zj 



JGI 



dzj 



dzi, 



(49) 



where {yi\j, ooj, zj) denotes the vector with I\J -coordinates from y, I -coordinates from 
z and other coordinates heing +oo. For instance, in case d = 2 we have 



oo pco o2 



yi -'Vi 



dzidz2 



(z) dz 



oo POO c)2 



yi •Jy2 



d^(P{z) 
dzidz2 



g{z) dz 



+ 



+ 



2/1 

oo 



?/2 



9{zi,y2) 



d^ 



dzi 



[zi,y2) -g[zi,oo) 



d^ 



dzi 



[zi,oo 



g{yi, z2)^—{yi, ^2) - g{oo, z2)-—{oo, Z2) 
0Z2 0Z2 



dzi 



dzo 



+ g{yi, 2/2)0(^1, 2/2) - g{oo, 2/2)0(00, 2/2) - g{yi, oo)0(yi, 00) + g{oo, 00)0(00, 00). (50) 

Assuming that for all y 



lim / z/(z, (iw) = 0, lim / i'{z,dw) = Q, (51) 

2— >— 00 / ^ 2—^00 / , 

J w>y J w<y 

equation (H71) rewrites in the equivalent conservative form 

{g{z) - g{y)) dz / — (;., dw) - / {g{z) - g{y)) dz / —{z, dw). 
-00 ^ w^y ^y J w<.y ^'^ 

(52) 



Proposition 2.6. Let L have form (H2l) with a hounded weakly continuous stochastic 
kernel v, so that L generates a C -Feller (i,e. its semigroup preserves continuous functions) 
jump process in R'^ that we denote X^ . Then L^^f^ is given hy (143 p . // the kernel v has 
continuous hounded mixed derivatives, so that 



dzi 



z, dw) 
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is again a hounded kernel (possibly signed) for any nonempty subset / G {1, ■ ■ ■ c?} (includ- 
ing the whole set {1, ■ ■ -d]), then L^^^^ can be rewritten as ( H6|) . Finally L^^^^ generates 
itself a C -Feller Markov process that we denote Y^ if and only if the following conditions 
hold: 

All mixed derivatives of orders from 1 to d — 1 of the jump rates are non-positive, i.e. 

^^{z,dw)<Q (53) 

for any proper subset I of {!,■ ■ ■ d}; and 

Qd 
' w>y 



/ ^ j—[z,dw)<Q, z>y. 



a'. <5*' 



// this is the case, the duality relation ( l20l) holds. 

Proof. Everything is proved apart from the criterion for the generation of a Markov pro- 
cess. To get it one only has to note that the operator J g{z)fi{y, dz) —(y{y)g{y) with given 
kernel n and function a is conditionally positive (and generates a process) if and only if 
the kernel /i is stochastic (i.e. positive), and that the kernels from various terms in ( H6|l 
are mutually singular, so that this positivity condition should be applied separately to 
each term. 

One completes the proof by the same argument as used at the end of the proof of 
Proposition 12.31 D 

Couple of remarks are in order here. Condition (!54l) is not very transparent. A simple 
particular case to have in mind is when the kernel z/ decomposes into a sum of kernels 
depending on all variables but for one, i.e. 

u{z, dw) = ^ z/j(2;i, ■ ■ ■ , Zj-i, Zj+i, ■■■ ,Zd, dw), 
j 

in which case the condition fl54|) becomes void (thus trivially satisfied). On the other 
hand, conditions f l53|) are easy to check. To visualize this condition it is instructive to 
observe that if g is a density of a positive measure on R"^, then the distribution function 



g(x) = / q{z)dz 

is positive, but has all mixed derivatives negative. Even more specifically, if u decomposes 
into a sum of kernels depending on one variable only, that is 

z/(z, dw) = y^ i^jizj, dw), 
j 

all conditions of Proposition 12 . 6 1 are reduced to an easy to check requirement that all rates 
/ i^j{zj, dw) are decreasing functions of 



Zj. 
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Let us note that the method of the calculation of dual used above can still be used 
for processes with a boundary. For instance, let us consider a process on R+ with the 
generator 

L(j){x) = I {(j){w) — (j){x))v{x,dw). (55) 

The operator V takes the form 



R4 



L'q{dz) = / q{x)i'{x,dz)dx — q{z)dz / v{z,dw) 
Jr+ Jn+ 



and the same calculations as above yield 



L^^^^g{y)= I g{z)dz —{z,dw)- g{z)dz I —{z,dw) 

Jo Jw>y '^Z Jy Jo<w<v ^Z 



- g{y) / z/(y, dw) + (7(0) / z/(0, dw), (56) 

that is, an additional term appears arising from additional boundary taken into account 
while integrating by parts. Under assumption (1511) . this rewrites in the equivalent conser- 
vative form 

ry r f),, 

L^'^^^giy) = / {g{z) - g{y)) dz / — (z, dw) 

Jo Jw>y C^ 

{g{z) - g{y)) dz f ^{z, dw) + f {g{0) - g{y))u{0, dw). (57) 

We assume strong smoothness condition for i/, which forces the dual Levy kernel to 
have a density. This is not necessary. Just assuming monotonicity of J^^ z/(z, dw) and 
/w< ^^(^5 '^^) (^^"^ thus the existence almost sure of non-negative derivatives of these 
functions of z), we obtain, instead of ( 158|1 . the formula 

/y r /'oo I* 

{g{z) - g{y)) 4 / i^{z, dw) - / {g{z) - g{y)) 4 / y{z, dw), (58) 
00 J w>y Jy Jw<y 

with similar modifications for (!57|) and analogously for (i- dimensional case. 

Let us mention the link with the theory of stochastic monotonicity. A Markov pro- 
cess Xf is called stochastically monotone with respect to Pareto ordering if the function 
P(Xf > y) is a monotone function of x for any y. Stochastic monotonicity is studied for 
various classes of processes, see [9], [10], [21], |23], [33], [32], [28] and references therein. If 
duality fl20l) holds, then Xf is obviously stochastically monotone, but, generally speaking, 
this condition is too weak to ensure duality, because stochastic monotonicity of a positive 
function on R'^ does not imply (apart from one-dimensional case) that it is the multi- 
dimensional distribution function for some positive measure. Therefore it is remarkable 
enough that for diffusion processes with generators 0351) the conditions of stochastic mono- 
tonicity and of the existence of Pareto dual coincide. Even for deterministic processes this 
is already not so, as for stochastic monotonicity of processes generated by operators (l29l) . 
bj are allowed to depend on other coordinates Xk (in a monotone way, see e.g. |9J and 
references therein to previous works). Stochastic monotonicity and related duality are 



well developed for Markov chains, see e.g. [2] and [3T], for birth and death processes and 
one-dimensional diffusions see [TT] . 

We assumed boundedness of all coefficients involved. This simplification leads to 
the most straightforward formulations that catch up the essence of duality. Of course, 
extensions to unbounded kernel rates, diffusion coefficients, etc, are possible under the 
conditions that ensure that all processes involved are well defined. 

2.6 Arbitrary Feller processes 

We have analyzed three classes of the generators L separately. But it is clear that if we 
consider a process with the generator being the sum of the generators of different classes, 
then applying conditions of the results above to each term separately will ensure that the 
dual to the sum is also conditionally positive and generates a process leading to the duality 
relation fl20|) . For simplicity, we shall give the corresponding result for one-dimensional 
Feller processes, but extension to higher dimensions is straightforward. For this case, the 
generators of the dual were obtained in [22j (which contains an annoying systematic typo 
with the wrong sign '— ' before the second term of flTTj) ) by approximating continuous state 
space generators by discrete Markov chains. The method of the present paper will give 
the same result without any technical restrictions used in [22] this yielding the complete 
characterization. 

Proposition 2.7. Let a Feller process Xf in Coo(R) have a generator 

Lg{x) = a{x)-—g{x) + b{x)—g{x)+ / {g{z)-g{x)-{z-x)g'{x)l\,^^\<i)u{x,dz) (59) 

with a,b & C^(R), a being non-negative, and with the weakly continuous Levy kernel v 
such that, for any y, conditions (ISTIl hold and the functions 



v{z,dw), — I v{z,dw) (60) 

u;>y J w<y 

are non- decreasing in z, for z < y and z > y respectively, so that their derivatives exist 
almost surely and are non-negative. Moreover 

'i-z<ydz iy{z,dw) + lz>ydz u{z,dw) (61) 

J w>y J w<y 

is a Levy kernel (it integrates min(l, {w — zY) and the integral 

ry+l 



lz<y{jy{y,dz)-^dz I u{z,dw)) + l^yy{u{y,dz) -d^ iy{z,dw)) 

w>y J w<y 



exists, at least in the sense of the main (or the Cauchy) value. Then the dual process Y^ 
exists (in the sense of fl20l) ) and has the generator 

d^ d 

Lg{y) = a{y)—g{y) + {a!{y) - b{y))-^g{y) 

+ / (^(^) - 9{y) -{z~ y)9\y)M^-y\<i)dz ( / y{z, dw) 

J —oo \J w>y 
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(giz) - g{y) -{z- y)g'{y)l\z-y\<i)d, ( / u{z, dw) 

y \Jw<y 



y+1 



+9'{y) / (z-y) 

Jy-l 



lz<y{i^{y,dz) + d^ I iy{z,dw)) + lzyy{iy{y,dz)-d^ u{z,dw)) 

w>y J w<y 

(62) 



Proof. Formula f l^ is obtained by combining fl55]) . f p7|) and fl?I]) . Conditions given ensure 
that the dual operator is well defined as a Levy-Khintchin type operator with variable 
coefficients. D 



Remark 7. ^45 shown in JEE/ and Theorem 5.9.2 of /^5]/ . conditions of stochastic mono- 
tonicity (monotonicity of functions fpl ) are sufficient for the operator (jHS]) to generate 
a Feller process, so that this condition can be dispensed with. 

As a corollary of Proposition 12.71 we can get now the full characterization of self - 
duality. 

Proposition 2.8. Let a Feller process X^ in Coo(R) have a generator ( l62l) . Then it is 
self dual (in the sense of (l20l) ) if and only if the following conditions holds: 

b{x) = a'{x)/2, dyu{y, dz) + dzU^z, dy) = 0. (63) 

In particular, if u has a density v{z.,w), which is differentiahle with respect to the first 
argument, then the second equation of (l63l) rewrites as 

^(y,z) + ^(z,y)=0. (64) 

Clearly, this condition is satisfied for v{y.,z) = g{\y — z\) with a smooth g, which corre- 
sponds to symmetric Levy generators. 

Proof. The condition on b follows from Proposition 12.41 The condition on z/ arises by the 
comparison of the integral terms of (1621) with (l62l) separately for y > z and y < z. D 

3 Stochastic /-duality from translation invariant / 

We have analyzed in some detail the duality arising from Pareto ordering. In general 
case explicit calculations are not always available. However, we shall propose here some 
general scheme for the analysis of translation-invariant /, that is / depending only on the 
difference of their arguments: 

f{x,y) = f{y-x), 

with some other function / that we still denote by / (with some ambiguity). 

Thus the operator F from ([5]) when applied to a measure Q with density q takes the 
form 

g{y) = {FQ){y) = [ f{y- x)q{dx), (65) 

JR'' 
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i.e. it becomes a convolution operator. It is then well known that under appropriate 
regularity assumptions, / is the fundamental solution of the pseudo-differential operator 
Lf with the symbol 

LfiP) = 7^, (66) 

where 



is the Fourier transform of /. 

Remark 8. In fact, by the definition of the fundamental solution, 



^'(7^)^W = ^W- 



which by taking the Fourier transform from both sides rewrites as 

Lf{p)f{p) = 1, 
as claimed. 

Hence g{y) from (!65l) solves the equation LfQ = q, so that F"^ = Lf. Of course, for an 
arbitrary /, the operator Lf can be quite awkward and the identification of the appropriate 
classes of functions q, g quite nontrivial. Let us consider the case when everything is well 
understood, namely the case of Lf being a Laplacian, or more generally, its fractional 
power. 

It is well known that the fundamental solution for the Laplace operator A in dimension 
c? > 3 is the function 

^^""^ = ~{d-2)a,^,\x\^-^' 

where a^-i is the area of the unit sphere in R*^. Hence the dual operator ([6]) takes the 
form 

rpDif) ^ ^-1 o T' o A, (67) 

and the generator for the corresponding dual semigroup becomes 

lD(f) ^ ^-1 o 2.' o A. (68) 

Let L be a diffusion operator of the special kind: 

Lg{x) = a{x)Ag{x) 
with a nonnegative bounded smooth function a{x). Then L' = Ao a{x) and thus 

lD{f) = A-^ o L' o A = L, (69) 

so that L is self /-dual. 

Noting that in dimensions d = 2 the fundamental solution for the Laplacian is known 
to be log |x|/27r we get the following. 
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Proposition 3.1. Let Xf he the Feller diffusion generated by the operator Lg[x) = 
a{x)Ag{x) in IV^ with a nonnegative bounded smooth function a{x). Then, for all x,y & 
R*^, we have 

Elog|Xf-y| = Elog|Xf-a;|, d = 2 (71) 

for c? > 3 and d = 2 respectively. 

Turning to the fractional Laplacian |A|"/^ in R*^ with a E (0,2), d > 2, let us recall 
that the inverse operator is given by the so-called Riesz potential 



|Ar-/2^(x) = I^gix) = -^ / 

Hd{a) Jr. 

r(«/2) 



gjy) dy 

Hd{a) Jn, \x-y\<'-'-' 
where 

^^(") = 2"^*'^r((.-„)/2)' 

see e.g. [IT]. Hence, the operator |A|"'^ is Lf for 

Let us consider a stable-like process generated by the operator 

Lg{x) = -a{x)\A\''/^g{x) 
with a positive smooth function a(x). Then L' = |A|"/^ o a{x) and thus 

L^(/) = |Ar"/2o^'o 1^1^/2^^^ (72) 

so that L is self /-dual. Thus we proved the following extension of Proposition 13.11 

Proposition 3.2. Let Xf be the stable-like process generated by the operator Lg{x) = 
a(x)\A\"^'^g(x) in R'^ with d >2, a G {0, 2] excluding the case d = a = 2 (for which ( ITTll 
holds), and with a nonnegative bounded smooth function a{x). Then, for all x,y E R*^, 

4 Stochastic duality for processes in R+ 

4.1 Reflected and absorbed diffusions in R_|_ 

We shall deduce some consequences from our general approach to processes on R+ that are 
dual in the sense ([3]). C^(R'^) will denote the space of k times differentiable functions on 
R*^ with all these derivatives vanishing at infinity. C^(R_|_) is the restriction of functions 
from C^(R) on R+ = {x > 0}. 

Consider a Feller process Xf on R generated by operator (I62p under the conditions of 
Proposition 12.71 assuming additionally that 
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(A) a G C^(R) and is an even function such that a{x) > 0, 6 G C^(R) and is an 
odd function (implying 6(0) = 0), the support of p is in R+ for x > and u{—x,dy) = 
Ri>{x,dy), where R denotes the reflection of the measure with respect to the origin (so 
that, by definition, f (l){y)Ri'{x,dy) = f (f>{—y)h'{x,dy)). 

Then, as is well known, see e.g. Theorem 6.8.1 in |23], the magnitude \X^\ is itself 
a Markov process on R+, also referred to as X^ reflected at the origin. Moreover, if the 
transition probabilities of X^ are pt{x, dy), then \X^\ has the transition density 

pT^{x, dy) = pt{x, dy) + Rptix, dy), 

and the semigroup T^'^ of |Xf | can be obtained from the semigroup Tt of Xf by the 
restriction to even functions. 

Remark 9. (i) Assuming that the kernel v is twice smooth would imply that the space 
C^(R) is an invariant core for X^ and consequently that the subspace of functions f from 
C^(R+) such that /'(O) = is an invariant core for \Xf\. (ii) If Xf were a diffusion, 
the process \Xf\ on R+ would he stochastically monotone by the coupling argument, see 
e.g. Sect 11,2 of 125]) and hence by Siegmund's theorem ^29] it had a Markov dual Y^ 
on R+ (in the sense ^) with absorbtion at the origin. In our case monotonicity follows 
from the construction of the dual below, which turns out to be given by a semigroup with 
a conditionally positive generator. 



Proposition 4.1. Under the conditions of Proposition 2.1, assumption (A) above and 
finally assuming that the measure z/(0, dw) is bounded, the dual process Y^ is a Feller on 
R_i_ absorbed at the origin and generated by the operator 

L^g{y) = a{y)^g{y) + {a'{y) - h{y))^g{y) + / (^(0) - g{y)v{(), dw) 
dy^ dy J^>y 

+ {9{z)- 9{y)-{z-y)g\y)l\,_y\<i)dA I u{z,dw) 

Jo \Jw>y 

{9{z) - g{y) -{z- y)g\y)l\,_y\<i)d, ( / z/(z, dw) 

y \Jw<y 

y+1 

lz<y{i^{y,dz) + d^ I u{z,dw)) + lz>y{u{y,dz) -d, I u{z,dw)) 

J w>y Jw<y 

(74) 



+9'{y) / (z-y) 

Jy-l 



The semigroup T/^ of Y^ is given explicitly by the formula 

/OO /"OO / /»oo C\ 

pT^{0,dz) + J^ g{x)i^J —prf{x,dz))dx. (75) 

Proof. Using ([8]) with F^^g{x) = —g\x) we get for g G C^(R_|. 



and hence 



{Tt9){y) = - dz g'{x)pr{x,dz)dx, (76) 

ly Jo 



/OO /"OO /"OO O 

pT^{0,dz) + J dxj g{x)^p7^{x,dz), (77) 
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yielding (!75|) as required. 

It is worth stressing that this formula implies the conservativity condition Tfl = 1 
(preservation of constants by Tf), because 

POO 

lim / pl^-^(x,dz) = 1 
by the Feller property and hence 

oo r\ / /"oo \ /"oo 

d^yj P7^^.dz)jdx = l-J p7^{0,dz). (78) 

Operators T/' form a semigroup by Proposition ll.2[ The form of the generator follows 
from (1571) . As it is conditionally positive, the semigroup Tf preserves positivity and 
preserves constants thus being a semigroup of a Markov process. Moreover, as also seen 
directly from (1751) . Tff{0) = /(O), so that the value at the origin is preserved meaning 
that this process is absorbing at the origin. D 

Remark 10. (i) Formula ( 1761) is valid only for g vanishing at infinity, and ( I75|) extends 
it (yields a minimal extension) to bounded functions on R+ . Plugging g = I into (175]) 
yields zero, not 1. (ii) The attempt to use integration in (1781) in opposite direction, at 
least when pt{x, dz) has a density pt{x, z), and using liuia-^oo Pi (^) 2;) = would give 

^^ \ I -^Pt'^x^ z)dx\ = - / pr^(0, z) dz, 

y \Jo (J^ / Jy 

which is different from the r.h.s. of (1781) . 

It is worth noting additionally that if a(0) 7^ and z/ = 0, then the subspace of 
functions g from C^(R+) such that 5'"(0) = is an invariant core for Y^ . In fact, the 
condition L^g{0) = (following from Tfg{{f) = g{0)) implies g"{0) = 0. On the other 
hand, if a(0) = and z/ = 0, then a{x) = ax^(l + o{l)),b{x) = bx{l + o(l)) as x — )■ 
with a > 0, 6 G R implying that is an unaccessible boundary point, so that Xf = |Xf | 
for a; > 0. In this case nothing comes out of the origin, so that p[^ (0,z) = for all 
z > implying that the first term on the r.h.s. of (I75|) vanishes and hence that is also 
unaccessible for Y^ (which follows also from its generator). In particular, if additionally 
b{x) = a'{x)/2, the process |Xf | is self-dual on R+. 

There is an extensive literature on the absorption - reflection link presented in Proposi- 
tion 14.1] mostly because of its natural interpretation in terms of ruin probabilities having 
important applications in insurance mathematics. For piecewise deterministic Markov 
processes it was obtained in [6] (see also [1]) and used effectively in |12j for assessing 
ruin probabilities via large deviations. Then it was extended to diffusions with jumps in 
[SO] , and to Levy processes in [5]. Our result is an extension of the corresponding result 
from [20], as we do it for arbitrary stochastically monotone processes. Our proof is quite 
different, as it is more elementary, using effectively only formula ([8]). 

4.2 Second dual and regularized dual 

Extension of the previous result to processes with a boundary from the right or with two 
boundaries is if course natural, see [S], but not quite straightforward. We shall clarify the 
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aspects of duality (even the definition has to be modified), needed for these cases reducing 
our attention to diffusions just for simphcity 

It is natural to ask whether the second dual coincides with the original process. For 
diffusions on R'^ this is in fact the case, as is seen from Proposition 12.31 However, for 
processes on R_|_ this dies not hold, as seen already from Levy's example of reflected 
Brownian motion. In fact, reflected BM cannot be dual to absorbing BM, as any dual 
process on R+ should be absorbing at the left end, that is at the origin, as seen directly 
from ([3]). However, the reflected BM is 'almost dual' to the absorbing BM in the sense 
that P(yi^ < x) = P(Xf > y) (with Y reflected and X absorbing BM) holds for all 
y ^ and all x. This suggests that the usual definition of duality imposes unnatural 
restrictions on the boundary. Consequently we shall give the following definition. Let X^ 
be a stochastically monotone process on [a, oo) such that P(Xf > y) is right continuous 
in X. A process Y^ on [a, cxd) will be called a regularized dual to a process Xf on [a, oo) 
if ([3]) holds for aA\ x > a,y > a, and the distribution ioi y = a is defined by continuity as 

P(F^<* <x) = PiYt"- <x) = lim P(F/ < x). (79) 

Remark 11. (i) One could also relax the condition for x = a defining P{Y^ < a) = 
lima,_^aP(y/ < x), but this would lead to the same result, as for usual definition, due to 
the right continuity of P(Xf > y) in x. (ii) If one only assumes monotonicity of the 
function P{X^ > y), it would become natural to define the dual distribution P{Y^ < x) 
as the right continuous modification of the function P{X^ > y). 

The following statement is now clear. 

Proposition 4.2. Under the assumptions of Proposition \4J\ the initial reflected process 
\Xf\ is a regularized dual to F/. Thus the second regularized dual to \Xf\ coincides with 

Remark 12. The usual (not regularized) dual of Y^ from Proposition \4-l\ is a rather 
pathological process Z^ , whose distributions coincides with that of |Xf | for z j^ 0, but the 
origin is an unattainable point without escape from it. Thus Z^ should be 'reflected from 
the origin' without touching it. 

Remark 13. Of course one can deal with reflected processes on R_ by introducing a 
symmetric notion of duality. Namely, for a process X^ on an interval o/R let us say that 
Y^ is its right dual, ifP{Y^ < x) = P{Xf > y) holds for all x,y (that is, it is the usual 
duality used above) and left dual ifP(Y^ < x) = P(Xf > y) holds for all x,y, which 
is equivalent to P{Y^ > a;) = P{Xf < y). Thus, by definition, Y^ is right dual to Xf 
if and only if X^ is left dual to Y^ . The theory of left dual processes on R_ (and their 
regularized version) is completely analogous to the theory of right dual process on R+. 
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